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ZJ ', Abstract 

^ ■ We consider the quantum radiation from a partially reflecting mov- 

ing mirror for the massless scalar field in 1+1 Minkowski space. Partial 
^ ■ reflectivity is achieved by localizing a 5-type potential at the mirror's 

H : 

5t , position. The radiated flux is exactly obtained for arbitrary motions 

as an integral functional of the mirror's past trajectory. Partial reflec- 
tivity corrections to the perfect mirror result are discussed. 



1 Introduction 

One of the most interesting prediction of the quantum field theory is the 
radiation produced by moving mirrors. A mirror can be practically seen 



as a barrier potential in the vicinity of a certain spatial boundary. Moving 
mirrors effects can be placed thus in the larger context of that produced in 
exterior classical fields. Essentially, the phenomenon expresses the fact that 
in presence of mirror-like boundaries in (non-uniform) motion real quanta 
will be excited out from vacuum fluctuations. 

For perfect mirrors, situation enjoys a certain technical definiteness. As 
infinite barriers are involved, the problem basically relies in assuring the 
quantum field to vanish on the corresponding boundaries. There are a lot of 
mathematical methods to deal with such problems, and a wealth of literature 
[0] following this line can be found. (Perfect reflectivity seems also to have 
been the preferred choice for discussing a series of most profound aspects of 
the phenomenon: the connection with the Hawking effect P]-P] and non- 
trivial situations implied by the possibility of generating negative energy 
fluxes Pj-lQ)- As concerns partial reflectivity, results are comparatively less 
numerous. Two main directions of research have to be mentioned here. One is 
that followed by Jaekel and Reynaud, based on the method developed in Refs. 
[^, |T^. Reflectivity properties of the mirror are described within a scattering 
approach, and the radiation is derived from the analysis of an S matrix 
depending on the mirror's trajectory. The method has a large generality, but 
it is associated with a major drawback: by its very nature it is ill suited for 
an accurate evaluation of the radiated flux density, as a local quantity, at 
arbitrary times. In addition, the analysis is carried out in frequency space, 
obscuring somewhat the trajectory dependence in coordinate space. Another 
direction is that of Barton et al. [jl^-|19[ (see also Ref. pO[). They proceed 



in the more conventional way of constructing the field operator via its modes 



expansion, tacking into account explicitly dielectric properties of the mirror. 

As a common feature, however, their results are systematically obtain within 

the framework of perturbation theory, in first order in the mirror's velocity. 

To our best knowledge, the only calculation allowing an exact result for 



the radiated flux is that of Frolov and Singh [£l| for spherical mirrors (in 
D-dimensional fiat space) expanding with a constant acceleration. Semi- 
transparency is realized by means of a 5-like potential. Our intention here to 
provide a two dimensional semitransparent mirror model based on a similar 
idea, allowing an exact treatment for arbitrary motions. 

We shall focus on the massless scalar field. The essential piece in our 
investigation is the construction of the (in) modes defining the field operator 
in Heisenberg image. We obtain them by introducing a set of space time de- 
pendent reflection and transmission coefficients into those of the Fulling and 



Davies model for perfect reflectors [^ . Matching conditions on the mirror's 



worldline force the coefficients to respect a first order differential equation, 
which determine them as integral functionals of the mirror's past trajectory. 
It turns out that calculations can be exactly carried out for the Wightman 
function and the renormalized energy-momentum tensor (we restrict here 
to the canonical expression). Qualitatively, one finds that the radiated fiux 
results as history dependent quantity, with past dependence extending in 
principle to minus infinity. Contributions of past motions appear exponen- 
tially damped with the proper time interval up to the emission moment, with 
the barrier strength as the proportionality factor. The perfect mirror result 
is recovered in the limit of infinite barrier strengths. 

The paper is organized as follows. In Sec. ^ we construct the quantum 



field. In Sec. ^ we obtain the in- vacuum Wightman function and tlie renor- 
Hialized energy-momentum tensor. A brief discussion of the near perfect 
mirror approximation (large barrier strengths or slowly varying motions) is 
also given. We end in Sec. ^ by drawing a parallel with extended charges in 
classical electrodynamics, relying on the history dependence of the emitted 
flux. Technical matters are referred to the Appendices. 

Natural units h = c = 1 are used throughout the paper. The metric 
tensor is gu = -Qxx = 1- 

2 The quantum field 

Let {t,x) denote the coordinates in 1-1-1 Minkowski space time and x{t) 
denote the trajectory of the mirror. We write the barrier potential as 

V{t,x) = ^^6{x-x{t)), 7=(l-(|)y', (1) 

where S is the Dirac distribution and a a positive constant corresponding to 
the barrier strength. The 7"^ contraction factor assures relativistic covari- 
ance of Eq. (|^) below. We proceed to construct the in modes for the field 
^{t,x) obeying 

(n + \/(t,x))<^(t,x) = 0. (2) 

We divide first the Minkowski plane into the left (L) and right {R) regions, 
corresponding to points {t,x < x{t)) and (t,x > x(t)), respectively. Let us 
write the in modes for the perfect mirror as 

Ut{u,v) = e-^-"^ - e-'-f^''\ U!^{u,v) = Q, (3) 



V^{u, v) = e-'^'' - e-^^3(^\ V^iu, v) = 0, cu > 0, (4) 

where L, R superscripts refer to the definition domains and null coordinates 
were used {u = t—x, v = t+x). Functions /, g are chosen such that identities 

f{v) = u, (5) 

g{u) = V, (6) 

are satisfied when {u, v) moves along the trajectory. 

Consider, for example, the U^ solution. e~'''^-^^'"' represents the totally 
refiected component of the incident right-moving wave e"*"^". Partial refiec- 
tivity requests introduction of a refiection factor, along with a transmitted 
component in R. We modify thus expressions (H) as 

U^{u,v) = e--"-7^;^(l;)e--^(^), (7) 

U^{u,v) = Tj'{u)e-'-\ (8) 

For V^'^ we similarly set 

Kf(w,i;) = e-^'^''-7^^(w)e-^'^^(^^ (9) 

V.^{u,v) = r^^{v)e-^-\ (10) 

The physical significance oiTZ^'^, T^'^ coefficients is clear. For points inside 
the L, R regions, the m, v dependence in Eqs. (|7|)-(0) automatically assures 

At the L-R boundary, for each pair of solutions (y^^, v!^)Lp=uy we impose the 
continuity condition 

^^ = ^t (12) 



and the matching of the first order derivatives, in conformity with Eqs. (|1|), 
(H). This is most conveniently obtained by considering the relativistic gen- 
erahzation of the stationary case x{t) =const. One finds 

e^.f/'^^V^ - e^.f/'^aVi + a y.^^ = 0, (13) 

where 6^,^ is the totally antisymmetric tensor in two dimensions {etx = 1) and 
f/'^ is the mirror's two velocity {U^U^ = !)• Equation (|12D implies 



7^^ + rJ^ = 7^f + r,^ = l, (i4) 

and Eq. (|l|) 

where all quantities refer to a given point on the trajectory. 

Let r denote the proper time of the mirror, and let us regard 7^, T 
coefficients as functions of r via u = u{t), v = v{t). Eliminating Tj^'^ in 
favour of TZ^'^ using Eq. (p!4D, one finds that Eqs. (|l^) and (|16D combined 
lead to the following differential equation for 7^^^'^ 

In^^'^r) + (-.c.Dj(r) + I) nt^ir) = ^, (17) 

where we introduced the Doppler factors (/? = dx/dt) 



Dtir) 



l±/3(r) 



(18) 



Plus and minus signs in D"^ correspond to R, L, respectively. For example, 
for uniform trajectories /3 =const the solution reads 

^^'^ = • n^\ „ ^ (19) 

—iujD^ + a/2 

having a transparent interpretation in terms of Doppler shifts as observed in 
the mirror's proper frame. 

The general solution can be written 

T^t'^'ir) = n^'^To) exp [zu;6^{r, tq) - a{T - To)) 

+ 1 £ dr' exp (tuSHr, r') - |(r - r')) , (20) 

with To fixed and 

S^{T„n)= rdr'D^ir'). (21) 

Jti 

With definitions (|18]) one actually finds for 6^ 



6+iT2,n) = v{T2)-v{n), (22) 

6-{T2,n) = u{T2)-u{n). (23) 

We let now tq -^ — oo. The first term in Eq. ( pO]) vanishesa, hence 

nt^ir) = "- £^ dr' exp (.^5±(r, r') -"-{r - r')) . (24) 

We take Eq. (^) as the definition for the reflection coefficients. Transmission 
factors follow from unitarity relation Eq.(0). This completes the derivation 
of the in modes. 



M-uLM. 



TZ^- {—oo)\ < 1 by the assumption of uniform velocity in the infinite past, see next. 



An observation is appropriate here: by the real part of the exponential 
factors, the coefficients are determined mainly by the motion in an interval 

< r - r' ~ a-\ (25) 

As an immediate consequence, consider a trajectory with constant velocity 
after some fixed proper time Tc. Then for sufficiently late times 

r - Tc > a'^ (26) 

non-uniformities in velocity before Tc can be ignored and one can set 

5Hry) = Df{r^){r-r% (27) 

which is equivalent to the stationary result ([19|). 

The Heisenberg field is constructed along the usual lines as 

V^(^) = /J j^ (a^ U^i^) + al V3{z) + H.c) , (28) 

for points z = (u, v) in the right region, and similarly with R —>■ L for z in 
the left region. The commutations for the creation-annihilation in operators 
are 

K, a^'1 = K, a^] = 2n{2uj)5{LO - c'), (29) 

and zero in rest. To assure that canonical commutations are satisfied, we 
shall suppose there exists a finite time to (which can be taken arbitrarily 
far in the past) such that the trajectory is uniform before to- Then, in the 
reference frame where the mirror is at rest at infinite past, the modes are 
given up to a certain moment by the (3 = coefficients (p!9|). Orthonormality 



and completeness can be easily checked in this case. It follows that canonical 
commutations are satisfied on some spacelike three surface, and thus, by Eq. 



they are respected in whole Minkowski space |24 



3 The radiation 

We are interested in the in- vacuum renormalized expectation values (T^;y)ren 
of the canonical energy-momentum operator 

Tf,y = d^^ifdyif - -g^^g"'^da^df3(f, (30) 

at points off the trajectory. We need the in-vacuum Wightman function 

D+{z,z') = {0,m\(P{z)(P{z')\0,m), (31) 

with the in vacuum defined by 

a^|0,in)=a^|0,in)=0, u > 0. (32) 

Let Dq{z, z') denote the free field Wightman function, and let's introduce 
the renormalized value 

D+ Jz, z') = D+iz, z') - D+{z, z'). (33) 



Then 



{Tuu{u,v))^^^ = ]im limdudu'D:^^^{u,v,u',v'), (34) 

{T^r{u,v))ren = l}m \im dydy> DXn{u,v, u' ,v') . (35) 



U'^U V'^rV 



Space time coordinates symmetrization proves unnecessary here, as Eqs. 
(|34D , (|35| ) will lead to real results. The mixed uv, vu components vanish 
identically. 

Let D^ stand for D^{z, z') when both z, z' belong to R, and analogously 
for D^. It is convenient to introduce Tu, Ty as the inverse functions of u, v 

Tu{u{t)) = T, r„(f (r)) = r. (36) 

Simple calculations yield (note that uj dependence in 7^^'^ makes the fre- 
quency integrations similar to that in the free field case) 

D+ = D+ + D+, + D+2, (37) 

where (e -^ 0+) 

(X fTu 

Dm = ^ dn In {{u{n) -u ~ ie){v{Ti) -v - it)) 



xexp-(ri -r„J 



+-— / " drx In ((u - uiji) - ie){v - v{ti) - ie)) 
Sir J-oo 



xexp-(ri -r„/), (38) 



and 



,2 



Dr2 = -TF~ / dTi / rfraln {{u{ti) - u{t2) - ie) 

iOTT J-oo J-oo 



x(f (ti) - v{t2) - ie)) exp -{n + T2 - r^ - r„'). 
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(39) 



The explicit form of D^ is of no relevance here. Tu is a shorthand for Tu{u), 
and similarly for r„/. D~j^ can be obtained from D^ with the substitutions 

(40) 



U ^^ V, m' <-> v' , Tu <-» Tjj. 



We shall take advantage in the following of the R — L formal symmetry above 
to refer to R quantities only. 

Renormalization amounts to ignore Dq in Eq. (^^. The double v deriva- 
tive (Bl) is trivial 



K{%v)), 



0. 



(41) 



The uu component ( pSf ) requires some effort. Calculations are outlined in 
Appendix A. One finds the v independent quantity 

1 l-P\ru 



{Tliu,v))r 



where 



with 



47r(l-/3(r„))^ 



jR ^ jR ^ jR^ 



XT^I 



r„ 



mru) = - 

4 J —OO J —OO 



'^u f^u 



dTidT2 






(42) 



(43) 



and 



T, 



Piiru) 



2 J- 



'^u r^i 



u I ' u 



xexp-(ri + T2 -2tu) 



dTidT2 dr^ dr2 —, ^ 7 7 



(44) 



xexp-(ri + T2 -2r„). 



(45) 
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The overdot represents differentiation with respect to proper time. Coinci- 
dence hmits in the integrands are discussed in Appendix B. 

Equations (|4^)-(p5D represent our main result. They exphcitly show that 



for finite barrier strengths the radiated energy-momentum density at a given 
point 

(r«(«,.;)),en = -(T,^(n,t;))„n = (T«(r„)),,„ (46) 



depends on the entire mirror's past historyn. Past dependence extends up to 
the retarded (emission) time r^, marking the intersection of the trajectory 
with the past hght cone at {u, v). As for TZ, T coefficients, exponential factors 
restrict the influence of past motions to an effective interval ~ a~^ before r^. 
Perfect reflectivity naturally results by letting a —>■ oo. To evaluate the 
limits, it is convenient to use 

lim a r dr'fir') expa(r' - r) = lim /(r'), (47) 

a^oo J —oo t'-*t 

t'<t 

which makes apparent disappearance of history dependence. For Eqs. (0), 
(^Sp situation is as follows. Let a denote the mirror's proper acceleration. 
Then coincidence limit (^) in T/^ can be rewritten (we assume the acceler- 
ation is continuously differentiable at r^, see Appendix B) 

lim Tf{Tu) = la{T^). (48) 

For TPj component, suffices to note that the double integral is multiplied by 
only one power of a, so that it brings no contribution. Combining Eqs. 



2f), one recovers the well known perfect mirror result [22, 23 



^A similar conclusion was reached, e.g., in Refs. [|^ |l5[ . 
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Suppose now that a is finite, but sufficiently large that no significative 
changes in velocity occur in the effective interval before r„. In other words, 
we are restricting to slow motions on a proper time scale ~ a~^. Then the 
brackets can be Taylor expanded around ti = r2 = r^, followed by a term by 
term integration. The result is the series 

T^ = To^ + E(-) T^^ (49) 

n=l ^'^^ 

with T^ coefficients completely determined as polynomials in a and its ffist 
n + 1 derivatives at Tu- T^ equals precisely the perfect mirror contribution 
(^Hf). Higher terms naturally describe semitransparency corrections. We 
give below the ffist two coefficients (for arbitrary orders, evaluation of the 
coincidence limits in the Taylor expansion proves quite laborious) 

T?- = -aa - -a, (50) 

6 6 

Tn = - — a a a aa -\ a. (51) 

^420 6 6 120 ^ ^ 

It is interesting to note that both corrections can be accommodated into the 
perfect mirror result by replacing a in ( ^8]) with an "effective acceleration" 
containing supplemental a, a dependent terms. 

We also point out the following aspect: one sees that corrections can 
generally both diminish or enhance the zeroth order flux. Consider, for ex- 
ample, a leftward accelerated trajectory on which the perfect mirror radiates 
a positive energy flux in R 

a < 0, d < 0. (52) 

Then if the flux is increasing or not decreasing too fast as seen in the mir- 
ror proper frame (i.e. a is not too large), the leading correction Eq. (pO|) 

13 



has always a diminishing effect. By contrast, for the absolute values of the 
negative energy flux radiated in L, the opposite is true. 

We end with a comment on non-radiative motions. It is clear that non- 
radiativity results as a history dependent property. Consider a fixed trajec- 
tory up to proper time tq. To find its non-radiative future extension, one 
should solve the integro-differential system 

r^(r) = r^(r) = 0, r>ro, (53) 

where t < tq trajectory plays the role of Cauchy data. Unfortunately, a 
general treatment is analytically rather difficult, if not impossible, given non- 
linearities implied. We shall contend ourselves to point out the special case 
when both integrands in Eqs. (pj), (^) vanish identically. Coincidence limit 
in Tf^ readily gives a necessary condition as 

a(r) = const. (54) 

One can easily check that this is also sufficient. Hence no radiation occurs 
for constant acceleration trajectories, and this is what one should expect by 
conformal invariance of the theoryH. An observation is in order here. It may 
be argued that uniformly accelerated trajectories do not respect the constant 
velocity condition for t -^ — oo as formulated in Sec. |^. From a practical 
point of view, however, this is not essential: the mirror can be considered 
stationary up to a certain moment, and only afterwards taken to accelerate. 
After a proper time large enough compared to a~^ the motion will effectively 
count as non-radiative. 



^A similar situation happens for D = 2 spherical mirrors in Ref. [gl[ 
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4 An analogy with extended classical charges 

We conclude by briefly discussing an analogy. History dependence of the 
emitted flux invites to a parallel with extended charges in classical electro- 
dynamics. Consider an arbitrarily moving charge with characteristic length 
r. Retardation implies that the electromagnetic flux at a certain point is 
determined by the motion in some finite length interval At ~ r. One is led 
thus to the correspondence 

r ~ a~^. (55) 

It is further relevant to turn toward a dynamical aspect. It is a common 
remark that the quantum backreaction force experienced by perfect mirrors in 
two dimensions |]22[| displays the same trajectory dependence as the radiation 
force p5| acting on pointlike charges. This makes the well known r —>■ 
pathologies, i.e. runaway trajectories and acausality, manifest too for perfect 
mirrors. 

Our observation concerning similarity expressed by (|55|) is as follows. It 
was shown in Refs. |]26|, ^7j that situation above changes when allowing for a 



transparency factor: more precisely, the unphysical behavior is absent when 
the transparency cut-off is sufficiently small compared to the mirror's massQ. 
Now, the cut-off in our model is some quantity of order a. These find a 
direct parallel in what happens for extended charges: the inverse of r plays 



^The result is based on the S matrix approach in Ref. iQ. There are certain hniitations 
in the proof, as approximations are made to hnearize the dynamical equation. In our 
model, linearization shows there are no unphysicalities if a/{AT:m) < 1, with m the mirror's 
mass. 
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an identical role, in the sense that unphysicalities disappear |^ provided the 
ratio between the electromagnetic self energy ~ r~^ and the physical mass is 
small enough. 
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Appendix A 

We sketch here the steps leading to Eqs. (|4^), p5|). Let 

ti,2{u, u') = dudu'D^-i^^iu, V, u, v'). (56) 

It turns out that ti, ^2 separately diverge for e -^ 0, u' —>■ u. An appropriate 
grouping of terms in their sum, however, makes divergences to cancel among 
themselves. We show below a way to do this. Observe that terms in ^2 can 
be organized in u, v dependent parts. We denote them as ^2; ^2; respectively. 
We further introduce a set of equivalent, but formally different expressions 

''2o ~ ''2/3 ~ ''27 ~ ''2? ''27 ~ ''2' {'^ ' ) 

which we define as follows. 
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With a series of integrations by parts with respect to ri, ^ can be cast 
into the form 



u _ __(^dTudTu^( r" , dui exp f (n - r„) 
" IGvr du du' \J-oo dri u{ti) — u' — ie 

-- r ^-1?^ r ^r,^"p'/";+"r."""""'n- (58) 

Similar integrations with respect to T2 lead to a second expression t2^. An 
integration by parts with respect to T2 in the double integral in Eq. (^) 
eliminates the single integration term. Let tl^ denote the resulting expression. 
Starting with t"^ and repeating the steps leading from ^2 to ^37) one obtains 

''27- 

We organize now the sum ti + 12 as 

K-t-2,) + it2a + t-2(^ + tl). (59) 

At this point one can set e = and let u' —>■ u, in each parenthesis the diver- 
gences getting cancelled in a manifest way. The two parentheses correspond, 
in order, to / and // contributions in {T^{u,v))ren- 



Appendix B 

We make here some observations concerning the coincidence limits in Eqs. 
(^11), (^51). Simple analysis shows that the quantities in the brackets are 



finite, provided the trajectory is at least of class C^, respectively C^. The 
coincidence limit in Tf^ is of direct relevance for the perfect reflectivity limit. 

For ri^2 ^ T it reads 

1 V (r) 1 u (r) 



6 v(t) 6 u{t) 
17 



(60) 



An equivalent way of writing T^jj is to eliminate the double derivatives 
dr^d-r^ by integrations by parts. Then coincidence limits involve only the first 
derivatives of m, v in T/*^, and the first two derivatives of u in T/^. The same 
applies for the resulting boundary terms at r^. The conclusion is that a C^ 
trajectory suffices to assure the continuity of the flux. Note that this is not 
the case for perfect mirrors. 

With the above choice, discontinuities in the acceleration manifest in 
T^i in the boundary term at Tu. Consider a trajectory with acceleration 
continuous everywhere excepting the moment Tq, where 

ii(ro+) - u{tqJ) = Auq ^ 0, ro± = tq + 0±. (61) 

One finds this entails the fiux discontinuity 

(rJJ,(w(ro+),t;)).en - (Tj^„(w(ro-),t;)).,„ = ^^^. (62) 

As a approaches infinity, the quantity above is progressively canceled by the 
integration terms in Tfj. In the limit, however, the discontinuity survives 
due to the Tf contribution (when it assumes a distributional form cf. (|60D). 

References 

[1] For a brief survey see W. G. Anderson and W. Israel, Phys. Rev. D 60, 
084003 (1999). 

[2] N. D. Birrell and P. C. Davies, Quantum Fields in Curved Space (Cam- 
bridge University Press, Cambridge, 1982). 

[3] P. C. Davies, J. Phys. A 8, 609 (1975). 



[4] W. R. Walker, Phys. Rev. D 31, 767 (1985). 

[5] R. Carlitz and S. Willey, Phys. Rev. D 36, 2327 (1987); ibid. 36, 2336 
(1987). 

[6] F. Wilczek in Black Holes, Membranes, Wormholes, and 
Super strings. Proceedings of the International Symposium, Woodlands, 
Texas, 1992, edited by S. Kalara and D. V. Nanopoulos (World Scien- 
tific, Singapore, 1993). 

[7] T. D. Chung and H. Verlinde, Nucl. Phys, B418, 305 (1994). 

[8] M. Hotta, M. Shino, M. Yoshimura, Progr. Theor. Phys. 91, 839 (1994). 

[9] P. C. Davies, Phys. Lett. 113B, 215 (1982). 

[10] D. Deutsch, A. C. Otewill and D. W. Sciama, Phys. Lett. 119B, 72 
(1982). 

[11] W. R. Walker, Class. Quant. Grav. 2, L37 (1985). 

[12] L. H. Ford and T. A. Roman, Phys. Rev. D 41, 3662 (1990). 

[13] W. G. Anderson, Phys. Rev. D 50, 4786 (1994). 

[14] L. H. Ford and T. A. Roman, Phys. Rev. D 60, 104018 (1999). 

[15] M. T. Jaekel and S. Reynaud, Quantum Optics 4, 39 (1992). 

[16] M. T. Jaekel and S. Reynaud, J. Phys. (France) I 2, 149 (1992); ibid. 3, 
1 (1993). 

19 



[17] G. Barton and C. Eberlein, Ann. Phys. 227, 222 (1993). 

[18] G. Barton, Ann. Phys. 245, 361 (1996). 

[19] G. Barton and C. A. North, Ann. Phys. 252, 72 (1996). 

[20] R. Gutig and C. Eberlein, J. Phys. A: Math. Gen. 31, 6819 (1998). 

[21] V. Frolov and D. Singh, Class. Quant. Grav., 16, 3693 (1999). 

[22] S. A. Fulling and P. C. Davies, Proc.R. Soc. London, A348, 393 (1976). 

[23] P. C. Davies and S. A. Fulling, Proc. R. Soc. London, A356, 237 (1977). 

[24] J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields (McGraw- 
Hill, New York, 1964). 

[25] F. Rohrlich, Classical Charged Particles (Addison- Wesley, Reading, 
MA, 1965). 

[26] M. T. Jaekel and S. Reynaud, Phys. Lett. A 167, 227 (1992). 

[27] M. T. Jaekel and S. Reynaud, J. Phys. (France) I 3, 339 (1993). 

[28] See H. M. Franca, G. C. Marques and A. J. da Silva, Nuovo Cim. 48A, 
65 (1978) and references therein. Proofs generally rely here too on the 
analysis of a linearized equation of motion. 



20 



